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Abstract 

>' 

cn i We present a calculation of 0{as) contributions to the process of t-cliannel single- 

top production and decay, which include virtual and real corrections arising from 
O ' interference of the production and decay subprocesses. The calculation is orga- 

t^^ ■ nized as a simultaneous expansion of the matrix elements in the couplings a^w, C(s 

^ . and the virtuahty of the intermediate top quark, {pf — mf)/mf ~ Tt/rrit, and 

extends earlier results beyond the narrow-width approximation. 

>'■ 
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1 Introduction 

Both the DO and CDF collaborations have recently announced the observation of single-top 
production at the Fermilab Tevatron at a significance of 5 standard deviations [llE]. This 
process represents a promising channel for the study of the charged- current weak interac- 
tions of the top quark, and will play a prominent role in the physics program at the LHC, 
where top quarks will be produced singly at large rates. Measurements of the single-top 
production cross section can be used to directly determine Vtb and to test the unitarity of 
the CKM matrix |3]. Furthermore, angular correlations of the products of the top-quark 
decay encode information on the spin structure of the Wtb vertex and on the production 
dynamics of the top quark |11|5] . The single-top production cross section also probes the 
bottom-quark distribution inside the proton, which at the moment is computed from light- 
parton densities rather than extracted from data. Therefore, this reaction represents a 
means of directly constraining the heavy-quark content of the proton. Finally, single-top 
production constitutes a background to a number of possible new-physics processes, most 
notably some channels important for Higgs boson searches. In view of all these consider- 
ations, a precise theoretical description of single-top production in hadronic collisions is 
highly desirable. 

In the Standard Model (SM) single-top production proceeds via three main hadronic 
channels, namely, t-channel production, qb — )■ q't or qb — )■ q't; s-channel production, qq' — )■ 
tb] and associated tW production, bg — ?■ W~t. At the Tevatron and at the LHC the t- 
channel production process has the largest cross section. In particular, at the LHC the 
s-channel production cross section and the tW production cross section are expected to 
be respectively 20 times and 3 times smaller than the t-channel process (see e.g. Ref. p]). 
Thus, in this paper we concentrate on the t-channel production mechanism. However, 
the distinction between t-channel and s-channel production is problematic beyond leading 
order and we will have to be more precise (in Section |2]) in describing how exactly we 
construct our observables and what we include in our calculation. 

Top-quark production, or the production of any unstable heavy particle, can be treated 
in several ways. The most straightforward way is to treat the top as a stable particle and 
ignore its decay. In this context, the cross section for the t-channel single-top production to 
next-to- leading order (NLO) in QCD was computed in Refs. [HE]. Later, a fully differential 
calculation was presented in Refs. |9l[T0]. NLO QCD corrections have also been computed 
for s-channel and associated tW production |9l|Tll[T2]. The full electroweak corrections in 
the SM and MSSM were computed in Ref. [13] for stable t-channel single-top production, 
and very recently for both t- and s-channel processes [H]. Finally, effects of soft-gluon 
corrections beyond NLO have also been studied [15]. 

Beyond the stable-top approximation, the one-loop corrections split in a gauge- invariant 
way into so-called factorizable and non-factorizable corrections [iGtlTT] . Factorizable correc- 
tions correspond to (on-shell) corrections to either the production or the decay part of the 
process. Thus, a way to simplify the calculation is to separately compute the corrections 
to the production and decay of an on-shell top. This approximation (sometimes referred 
to as the improved narrow-width approximation) allows the inclusion of realistic cuts on 



the decay products of the top. NLO QCD analyses in this framework for the semileptonic 
top decay were pubhshed in Refs. [T8142T] . 

To our knowledge, none of the presently published works on t-channel single-top pro- 
duction include the effects of interference between real radiation in production and decay 
or virtual corrections connecting the two subprocesses. A study of these non-factorizable 
contributions has been presented for s-channel single-top production [22], and for tt pro- 
duction [23]. These corrections are known to be very small, for observables which are 
inclusive enough in the invariant mass of the top quark pll|25] , due to large cancellations 
between virtual and real contributions. However, there is, a priori, no reason why this 
should hold true for arbitrary observables, especially if they involve kinematical cuts that 
could, in principle, spoil the delicate cancellation of real and virtual contributions. 

In this paper we want to assess the effect of these production-decay interferences in 
t-channel single-top production at NLO in the QCD coupling constant, a^. Hence, we 
are interested in resonant top-quark production with pf ^ mf, but pf ~ mf or, more 
precisely, p"^ — m^ ~ fntTt, where Fj ~ 1.4 GeV is the width of the top. While the effect 
of these "off-shell" corrections is expected to be very modest for the total cross section, 
we are particularly interested in distributions that are related to the measurement of the 
top-quark mass, rrit. For a reliable mass determination with an error 6mt < Ft, the non- 
factorizable corrections have to be under control. In this respect, we also view the current 
calculation as a preparation to apply our method to top-quark pair production. 

We will neglect quark-mixing effects and treat the bottom quark as massless throughout, 
using the 5-flavour scheme. The importance of bottom-quark mass corrections and the 
relation between the 4- flavour and 5-flavour scheme has been studied in Refs. [26|[27]. 
Furthermore, our calculation does not include any effects due to parton showers. The 
matching of the NLO QCD result with parton shower Monte Carlo was implemented in 
MC@NLO [2H] and in POWHEG [21]. 

The calculation is organized as an expansion in the virtuality of the top quark, p^ — mf, 
in a way similar to the pole approximation [301 [31], and employs effective-theory (ET) 
inspired techniques analogous to the ones used in Refs. [3214M] . These are based on splitting 
contributions to the amplitude into so-called hard and soft parts using the method of 
regions [5^15^ , thereby extending the separation between factorizable and non-factorizable 
corrections beyond NLO. The hard part can be identified with the factorizable corrections, 
whereas the soft part corresponds to the non-factorizable contributions [37]. This approach 
has the advantage of providing a gauge-invariant resummation of top-quark finite-width 
effects. Furthermore, the expansion in pf — mf allows for an identification of the terms 
relevant to the achievement of a given target accuracy before the actual computation, 
leading to a significant simplification of the calculation. The method has been discussed 
in detail for a toy model [33] and can easily be adapted to our case for the tree-level and 
virtual contributions. For the real corrections this is more problematic and we will not be 
able to follow a strict ET approach in this case, but will have to combine ideas from the 
effective theory with a standard fixed-order approach. As will be shown in this paper, this 
results in straightforward calculations of the contributions that are expected to be relevant 
for phenomenological applications at hadron colliders. 



The outline of the paper is as follows: we start with a general description of our method 
to deal with resonant particles at hadron colliders in Section |2l While we will concentrate 
on the process at hand, the discussion is meant to be general enough to be easily adapted 
to other processes. We will also be more precise in describing the observables we are 
interested in and the accuracy we are aiming for. In Section [3] we first give explicit results 
for the amplitudes needed and details of the computation. We then discuss a series of 
successive approximations to the exact cross section, which relate to previous results in 
the literature, and illustrate the cancellation of infrared singularities in the various cases. 
Numerical results for the cross sections and distributions will be presented in Section HI 
Finally, in Section [5] we summarize and give an outlook on further possible applications of 
our method. 

2 Method 

2.1 Setup of the calculation 

In order to include production-decay interference effects, the narrow-width approximation 
has to be relaxed. In particular, the top quark cannot be treated as a stable particle. 
Physical observables must be computed for final states containing only long-enough lived 
particles, which in the case of interest in this paper are represented by the products of the 
top-quark decay. Considering the LHC, we are interested in the process 

p{Pi)p{P2) ^ UPb)W+{pw)X, (2.1) 

where Jh is a 6-quark jet and X stands for an arbitrary number of further jets, as long 
as they do not originate from a 6 or 6 quark. Rejecting 6-quark jets (naively) excludes 
contributions from the s-channel process. Furthermore, we do not allow a second W in the 
final state in order to exclude associated production, and we insist on a positively charged 
lepton (from the W decay) to exclude single-t production. Some of these constraints are, of 
course, questionable from an experimental point of view and most could easily be avoided. 
But this is a minimal, more or less realistic setup that allows us to discuss the inclusion of 
non-factorizable corrections. 

The most important constraint we make on the final state is that the invariant mass of 
the iy"*"Jb pair is close to the top-quark mass, i.e. that the top is resonant. More precisely, 
we requirqj {ph + PwY ~ "^t ~ "^tfi ~ o^-ew'^'i ^ '^t • ^^ is well known, in this case 
strict fixed-order perturbation theory breaks down due to the kinematic enhancement of 
formally higher-order corrections. These corrections have to be resummed in a consistent 
way. 



^If taken at face value this would correspond to an invariant-mass window Apt ^ Tt around the top- 
quark mass. In fact, a numerical study reveals that the suppression of non-resonant configurations is 
already effective for much looser invariant-mass cuts. 






Figure 1: Tree-level Feynman diagrams for the process qb — )• q'bW. The figure shows 
both purely EW contributions, diagrams (a) and (c), and mixed QCD-EW contributions, 
diagrams (b). The semileptonic (on-shell) decay of the W is understood. 



To illustrate this, consider the process (12. ip at partonic tree level, where we have to 
compute 

qiPi) KP2) -> q'iPs) Kp^) W+ipw) -> q'iPs) h{p^) l+{p^) vi{p^) , (2.2) 

with the initial parton, g, being a light quark (m, c) or antiquark (d, s). Accordingly, q' 
is either a quark ((i, s) or antiquark (u, c) respectively. For the purpose of the current 
discussion the h quark can be identified with a b jet, J^, with momentum pb = p4. The 
decay of the W~^ — t- l^ui is described in the improved narrow-width approximation. This 
considerably simplifies the calculation, and yet allows for a non-inclusive treatment of the 
leptons coming from the VT-boson decay. 

The matrix element for the partonic process (12.21) can be computed from the Feynman 
diagrams shown in Figure [TJ These can be divided into two classes: resonant diagrams 
(diagram (a)), in which the Wb pair originates from the decay of an internal top quark; and 
non-resonant or background diagrams (diagrams (b) and (c)), that do not contain inter- 
mediate top-quark lines. The latter are subdivided into electroweak-mediated (diagrams 
(c)) and mixed QCD-EW diagrams (diagrams (b)). It is important to note that only the 
sum of all electroweak diagrams (diagrams (a) and (c)) is gauge invariant, though, strictly 
speaking, only diagram (a) describes the production of a single top quark. Obviously, the 
sum of all QCD-EW diagrams is separately gauge independent. 



The fixed-order tree- level amplitude, A , can be written as 
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where /C(pi) and J{pi) are functions of the external momenta, pj, i = 1,...,6. The first 
term in Eq. ( 12.3^ describes resonant contributions, whereas JT accounts for non-resonant 
diagrams. Eq. (12. 3p has a pole at p^ = mf that is related to the breakdown of fixed-order 
perturbation theory mentioned above. When an intermediate virtual top quark approaches 
the mass shell, a subset of formally subleading corrections to the top-quark propagator are 
enhanced, and must be resummed to all orders: 



+ mt] 



-)■ 



rrit 



Pt 



mi 



Pt 



mt 



E 

n=0 



-i^ti^t) 



i(^t + mt] 



Pt 



mi 



(2.4) 



f^ti'^t) denotes the sum of (renormalized) one-particle irreducible corrections to the top- 
quark two-point function, and contains an imaginary part of order ttetu^^t that, upon re- 
summation, regularizes the propagator. However, mixing different orders in perturbation 
theory leads, in general, to violation of gauge invariance and unitarity, which are guar- 
anteed only for strictly fixed-order calculations and for the full amplitude. Therefore, a 
meaningful gauge-invariant expansion of the matrix element in p^ — m,f requires addressing 
the issue of resummation of finite-width effects. 

In recent years, several frameworks for the gauge-invariant calculation of production 
cross sections for unstable particles have been proposed, from ad-hoc procedures to treat 
weak gauge bosons [38| , to more general methods based on complex renormalization [39tH0] 
or modified perturbative expansion in terms of distributions [UHlS]. Here we adopt an 



effective-theory inspired approach [33] which, in this case, can be considered as an exten- 
sion of the pole- approximation scheme j5Ull^ . As we will see, this provides a consistent 
framework for a systematic gauge-invariant expansion of the amplitude around the com- 
plex pole of the full top propagator at /i^ = m^ — imtVt. Note that we have used the same 



notation, m^, for the mass parameter in Eq. (12.41) and for the real part of the complex 
pole. Strictly speaking, the first one is the scheme-dependent renormalized mass, mt^r, and 
the latter the physical pole masci, and in a generic renormalization scheme they differ by 
an amount proportional to a^, as shown in Section 1X51 However, in the on-shell scheme 
adopted in this paper, the two masses coincide up to corrections which are beyond our 
target accuracy. Therefore, in the following we refrain from introducing more than one 
mass parameter, except for Section 13.31 where it will be necessary to explicitly distinguish 
the renormalized mass, m^ r, from the pole mass, ntt. 

For the tree- level amphtude, Eq. (12.31) . the pole expansion reads 
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^It is well known that for a mass determination with an error Srrit < Aqcd, a mass definition other than 
the pole mass has to be used [44ll45] . However, in this paper we ignore effects of the order Aqcd ^ nitS. 



where we have introduced A^ = pf — fif, and 1 + SRt denotes the residue of the full 
propagator at pf = fif. The leading resonant contribution is encoded in ]C{pi;pf = fif). In 
dIC/dpf the resonant propagator has been cancelled by a A^ arising from the expansion 
of }C{pi) around p^ = /i^, while Jipupj = [if) represents the leading contribution of the 
truly non-resonant diagrams. The ellipsis denotes higher-order contributions suppressed 
by additional powers of A^. These can, in principle, be computed to any order, but in the 
following discussion only the terms shown in Eq. fl2.5p will be relevant. 

Beyond the leading approximation, A}'^'^'^ ~ IC{pi;pf = fif)/At, the consistency of the 
expansion fl2.5p requires the inclusion of loop-corrections. The kinematic expansion in the 
parameter, Aj, must be combined with a standard expansion in the coupling constants, 
as = g1/{ATi) and aew = olwli^^)- Fo^ counting purposes, in the following we refer to the 
three expansions parameters collectively as 5 and assume the relative scaling 

(5 ~ a^ ~ Oe^ ~ ^ . (2.6) 

mi 

We are thus led to write the tree-level amplitude for the process (12. 2 p as 

^t- = (531(542 [gL Af^) + qL ^Sof ^ + •••)+ T^in2 de^gl A(^'2) . (2.7) 

The superscripts denote the order of the couplings which multiply the amplitude, while the 
subscripts denote the order to which the propagator, A(, appears within the amplitude, 
i.e. ^[n^'" has a prefactor g^g^ with Am'" ~ Aj. The coupling from the decay of the 
W is not counted. A missing subscript indicates that the amplitude does not contain 
a Ai propagator, i.e. A^"^-'^) = A^^f\ Thus, ^L ^(-?) ~ ^^''^^ gey^gl A^^^^~^ ~ 5 and 

dlw ^(0) ~ ^^ • Terms suppressed beyond 5'^''^ are indicated by the ellipsis. In the notation 

of Eq. (12. 5p . (yfg^y4;J-^; corresponds to }C{pi;pf = fif)/At, while g^^A^^^. includes both 

the dlC/dpi and J terms. Clearly, gewgl A*^^'^^ receives contributions from non-resonant 
diagrams only. 

Note that the QCD-EW contribution, A^^''^\ is usually considered to be a background 
to single-top production. However, the final state is identical and, therefore, A^^''^'> has to 
be included in Ji'^^^. In principle, there could be interferences between A,^A and A^^''^\ 
however, due to colour they vanish at tree level. For the colour-averaged square(l amplitude, 
M^""" = ]^ Ec l^*'''*'^, we obtain 
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where, as usual, Nc = ?> and Cp = 4/3. The (first) leading term of Eq. (12. 8p scales as 5, 
whereas the other two terms scale as 5^ and represent a correction of order 5 ~ 1% to the 
leading contribution. All other terms are further suppressed. 

Our aim is to compute all contributions to the cross section up to 0{5^''^). According 
to our counting, Eq. (12. 6p . this requires, beside the leading tree-level amplitude, the calcu- 
lation of A,^-^.^ the 0{as) one-loop corrections to the leading resonant contribution. For 



the squared amplitude this leads to a contribution 

M^^°=^1^.^2Re(Af?![Afj?J]*), (2.9) 

which we will refer to as NLO. Along with the virtual corrections we also have to include 
real corrections, g6 — )■ q'hW^ g. There are also gluon-initiated processes at NLO, namely 
gh — )■ q' h W^ q and qg — )■ q' h W^ h. The former is fully included in our calculation, 
whereas the latter deserves special discussion as it mixes t-channel and s-channel single- 
top production. Writing the amplitude for g(pi) g{j)2) -^ l'ij>z) KP4:) W^{Pw) b^pj) as 

<r = 9s qL {n^hi Al + Tllb,, Aji) , (2.10) 

we first note that upon squaring J^IT there are no interference terms, due to colour. 
Conventionally, the term ~ l^flP is included in t-channel single-top production, whereas 
the term ~ l^ggP is considered to be a higher-order correction to s-channel single-top 
production. We will follow this convention but stress that a fully satisfactory solution 
requires the simultaneous inclusion of both processes, which we reserve for future work. 
In summary, we need to compute the tree-level amplitudes for 

uh^dhW"^ db^ubW^ (2.11) 

ub-^dbW+g db ^ ubW^ g (2.12) 

ug^dbW+b dg^ubW+b (2.13) 

gb^dbW+u (2.14) 

and the QCD one-loop corrections to the leading (in At/rrit) part of process (12. lip . Of 
course, we will also have to include the processes with {u,d} — >■ {c,s}. Note that the 
various processes containing a gluon are related by crossing. We stress once more that 
the semileptonic decay, W~^ — )■ i^ue, is taken into account in the improved narrow- width 
approximation. 

2.2 Loop corrections in the effective-theory approach 

Detailed discussions and applications of effective-theory methods to the calculation of 
processes involving unstable particles can be found in the literature [32l|33l|371ll6]. Here 
we will restrict ourselves to a discussion of the points directly relevant to our calculation. 
The main idea of the ET approach is to systematically exploit the hierarchy of scales, 
[p^ — m^) ~ rritTt <^ m^, by integrating out degrees of freedom with virtuality ~ mf. 
In doing this, the Lagrangian of the underlying theory is rewritten as a series of gauge- 
invariant operators multiplied by matching coefficients, which are determined such that 
the ET reproduces the results of the underlying theory up to a certain approximation. 
The matching coefficients are guaranteed to be gauge independent and they contain the 
information on the degrees of freedom that have been integrated out. 



Once the hard part is integrated out, we are left with degrees of freedom with virtuahty 
much smaller than mf. In general, we have to take into account several different such 
degrees of freedom, but in the case at hand it is sufficient (to the accuracy we are aiming 
at) to consider only soft modes. 

Within this picture, the leading contribution to the process (12 .2^ is to be viewed as the 
production of an on-shell top, the propagation of a soft top, and the subsequent decay of an 
on-shell top. This is to say that the matching of the full theory onto the ET is done on-shell. 
It is important to note that on-shell in this context means pf = /x^, i.e. on the complex 
pole of the two-point function. Thus, the matching coefficients are generally complex. The 
propagator of the soft top is derived from the bilinear operator in the effective theory and 
includes the resummation of the self-energy insertions. 

There are several sources of higher-order corrections. First, there are 0{as) and 0{aew) 
corrections to the matching coefficients. These correspond to loop corrections to the pro- 
duction and decay of a stable top and have been taken into account in previous calculations. 
Furthermore, there are subleading corrections to the bilinear operator, i.e. corrections to 
the propagator. There are also subleading operators in the ET that do not refer to a 
top quark at all. These operators reproduce the effect of background diagrams, as well 
as subleading effects of resonant diagrams, and correspond to the terms J^ and dK/dpf 
in Eq. (12. 5p respectively. Finally, there are the contributions from loop diagrams in the 
effective theory, i.e. loop diagrams with the still dynamical soft degrees of freedom. These 
contributions correspond to the non-factorizable corrections and link the production and 
decay parts of the process. 

From a technical point of view, the calculation of the aforementioned corrections is 
achieved by using the method of regions [35l|36] and computing the hard and soft parts of 
loop integrals. The hard part is defined by expanding the integrand under the assumption 
that the loop momentum, q^, scales as q ^ nit. In what follows we will always suppress 
powers of mt in the scaling relations and write the scaling of the hard modes as g^ ~ 1. 
The soft part is obtained by expanding the integrand of a loop integral assuming that the 
loop momentum scales as g ~ (5. 

A strict application of the effective theory would require the introduction of collinear 
fields for the massless external fermions and a heavy-quark field for the top quark, as well 
as soft and collinear gauge-boson fields. However, we refrain from doing this because for 
the real corrections we will have to deviate from a strict ET approach anyway. We will, 
however, make extensive use of the gauge-invariant separation of the one-loop contribution 
into hard and soft parts, and the associated counting rules to obtain a result that is gauge 
invariant and reproduces all terms to the desired accuracy with a minimal amount of 
computation. 

The one-loop diagrams necessary to achieving (9(5^/^) accuracy are shown in Figure [2l 
Apart from the top self-energy insertion (diagram (EJd)), at the accuracy we are interested 
in we have to deal with QCD one-loop diagrams only. Due to colour, the one-loop diagrams 
containing gluon exchange between the upper and lower quark lines result in a vanishing 
contribution to the amplitude squared and are therefore not shown in the figure. A detailed 
discussion of the computation of the various diagrams will be given in Section 13.21 Here 
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Figure 2: Virtual QCD corrections to t-channel single-top production at leading order in 

Ai/mt. 



we briefly mention which parts of the diagrams in Figure |2] enter our calculation. 

As explained in detail in Section I3l2l taking the leading hard and soft parts of the QCD 
self-energy insertion (diagram (2a)) we obtain the scalings gl^ag 5^"^ ~ 1 and glyjds 5~^ ~ 5 
respectively. For the soft part this corresponds precisely to a correction we are interested 
in. The hard part, however, seems to be superleading, i.e. the one-loop correction is 
enhanced compared to the tree-level amplitude (which scales as glyj/5 ~ 5^/^). However, 
as we will see below, in the on-shell scheme used in this paper this superleading contribution 
is cancelled precisely by the counterterm. In a generic renormalization scheme this is not 
the case, and the leading hard part of the self-energy has to be resummed and enters the 
deflnition of the complex pole, /x^, as explained in detail in Section [3731 The same applies to 
the two-loop QCD self-energy. An insertion of the leading hard part results in a diagram 
scaling as gl^otl S~'^ ~ S^^'^ and, in general, has to be resummed. The subleading hard 
terms are suppressed by at least one factor p'f — mf ~ 6 compared to the leading hard 
contribution and, therefore, result in contributions g^w^^s ^^^ ~ S, i-e. of the same order as 
the leading soft contribution. Thus, the subleading hard part of the self-energy diagram 
has to be included (but not necessarily resummed) independent of the renormalization 
scheme adopted. 

To obtain the scaling of the EW one-loop self-energy insertion we simply have to replace 
tts by ttew Thus, diagram (2b) scales as 5^^"^ and 5^^'^ with hard and soft EW self-energy 
insertions respectively. Therefore, the hard part is not suppressed with respect to the tree- 
level amplitude and has to be resummed. The soft part contributes beyond the accuracy 
of our calculation and can be neglected. Obviously, the resummation of the hard part of 



the EW self-energy insertion corresponds to the resummation indicated in Eq. (12. 4p . 

For diagram (2c), the decomposition into hard and soft part is trivial in that the soft 
part results in scaleless integrals and therefore vanishes, whereas the hard part corresponds 
to the full diagram. This is not surprising as the loop correction in this diagram is not 
affected by the instability of the top quark. 

Diagrams (2d) and (2e) are more interesting and both behave in a similar way. The soft 
and hard parts scale as g^y^Us ^~^ ~ <^ and contribute at NLO. From an ET point of view, 
the hard part of diagram (2d) contributes to the matching coefficient of the production 
operator (or to the decay operator in the case of diagram (2e)), whereas the soft part is 
reproduced by a loop diagram in the effective theory. 

Finally, performing a similar expansion for the box diagram (2f), it can be seen that 
the hard part scales as 5"^ and thus contributes beyond the accuracy of our calculation. 
The soft part, however, scales as 5 and must be included at NLO. This illustrates the 
simplifications that can be achieved in the calculation. Rather than having to compute 
a full box diagram with several scales, we end up computing only the soft part, which is 
much simpler. 

2.3 Real corrections 

For the tree-level and virtual calculations, an ET approach can be used in a straightfor- 
ward manner. If one is only interested in the total cross section, real corrections can also 
be tackled in this way by relating the total cross section to the imaginary part of the 
forward-scattering amplitude. However, here we are interested in computing an arbitrary 
infrared-safe observable. Hence, we want to compute the real corrections without explicitly 
specifying the observable. 

It is clear that, to a certain extent, this is in confiict with an ET approach. An effective 
theory relies on making all scales explicit. Since the definition of the observable itself 
can introduce additional scales, it is not possible to follow a strict ET approach. In this 
subsection we describe how we deal with the real corrections, making sure that we keep all 
terms to the desired accuracy, for a general observable. The only assumption we make on 
the observable is that it does not introduce another hierarchy of scales. We stress that this 
assumption has to be made for any fixed-order calculation because a large ratio of scales 
usually introduces large logarithms which have to be resummed. 

To begin with, we need the real matrix element squared for the process (I2.12p . We show 
the diagrams containing a top-quark propagator in Figure O There are additional diagrams 
corresponding to diagrams (b) and (c) of Figure [T] with an additional gluon radiated off, 
however, they are suppressed and contribute only beyond the accuracy we are aiming at. 
The real corrections have an additional factor of a^, which allows us to work at leading 
order in A^/mt. 

In the case of the real corrections it is not, a priori, clear what the correct expansion 
parameter is. Diagrams (a)-(c) have a resonant top propagator for {pw + p^)"^ ~ m^, 
diagram (e) has a resonant top propagator for {pw + Vi^ VgY — "^t; and diagram (d) is 
resonant in both kinematic configurations. Depending on how the final state partons are 
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Figure 3: Real QCD corrections to resonant t-channel single-top production at leading order 
in /^t/mt. 



combined to jets, both regions can be relevant and have to be taken into account. This 
can be achieved by a shght modification of the usual subtraction method. 

In order to isolate infrared singularities, real corrections are usually computed by sub- 
tracting from the real, n + 1 parton matrix element squared, M^+i, a term M^'^f^^. that 
approximates the full matrix element in all singular regions. This term has essentially 
n-parton kinematics. It is added back and a partial phase-space integration is performed 
analytically to recover the infrared 1/e poles explicitly. The kinematic configurations asso- 
ciated with M^"^-^x correspond either to a gluon being soft or two partons being collinear. 

Thus, for M^',"f-^s it is always clear whether the gluon is combined with the 6-quark or not 
and, correspondingly, what the appropriate expansion parameter is. Therefore, we can 
expand the term that is added back and write for the real corrections 



/ 



(i$n+lM 



n+l^Wn+1 



d^ 



n+l 



K'L 



n+l 



Kr^)) + 



d<D„,i M!^°^. 



^n+l 



n(+l) 



(2.15) 



rf$ 



n+l 



M 



n+l 



n{+l) 



+ 



a^n+l ^^2n(+l) 



The explicit infrared poles in M,f^^. now match the infrared poles of the virtual cor- 
rections. The error introduced in the second line of Eq. fl2.15p is suppressed by a factor 
as 5 relative to the Born term and, therefore, is beyond the accuracy we are aiming at. It 
corresponds to the error introduced in expanding the virtual corrections. Treating the real 
corrections in this way, we ensure that the expansion of the real and virtual corrections 
in At/nit is performed in a consistent way, such that the infrared poles always match. 
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We have implemented this procedure in two independent programs, using the FKS |17] 
and dipole subtraction [48j methods. The agreement of the results produced with the two 
programs provides us with a useful check of the implementation. 

3 Helicity amplitudes for single-top production 

In this section we present all the helicity amplitudes necessary for the computation of the 
0{S) and (9(5^/^) contributions to the cross section of the process (12. ip . These consist of 
the leading resonant tree-level diagram (la) shown in Figure [Tj which scales as 5^/^, and of 
virtual and real QCD corrections to this diagram (including the gluon-initiated processes 
Eqs. (I2.14P and (I2.13P ). which scale parametrically as S^^'^as ~ 5. The resummation 
of finite- width effects, illustrated in Eq. (12. 5p . also requires the calculation of one- loop 
electroweak self-energies, as discussed in Sections 13.21 and 13.31 Any other contribution 
is suppressed by at least 6 compared to the leading tree-level process and is beyond the 
accuracy pursued in this work. The inclusion of real corrections in our formalism has 
been discussed in the previous section, and here we limit ourselves to giving a list of 
the relevant amplitudes. For the virtual corrections we will explain, in some detail, the 
application of the method of regions to the computation of the expansion in 6 of loop 
integrals. Throughout this section, and the rest of the paper, we adopt the helicity notation 
introduced in Ref. [19],[50] and make use of the following abbreviations: Sij = {pi +pj)^. 
Pi = P4 + P5 + P6, Dt=Pt - ml and A^ = pj - /if. 

3.1 Tree-level amplitude 

As pointed out in Section [2711 of the tree-level diagrams shown in Figure[T]only the resonant 
one, diagram ([Iti), contributes to the amplitude at C((5^/^). Given its purely electroweak 
nature, the only non-vanishing amplitude is the one for the helicity configuration uibL — )■ 
dhbh^^T^L and it reads 

.3 4(3,0) _ 4 r^^[12](46)(3|4 + 6|5] 



^'^'""^(-^^ ^^"'VMv.Fh. (3i3 + M2.)A, • ^^-^^ 



In Eq. (13. ip . the extra power of the coupling constant, Qew, and the prefactor ^J^: /MwVy/ 
arise from the inclusion of the W^-boson decay in the improved narrow-width approxima- 
tion; l/((p^ - M^f + Ml^Vlr) ~ 'k/{MwTw)KPw - ^w)- Note, however, that this 
prefactor does not modify the scaling of the amplitude since Vy/ oc gl^. The complex pole, 
/if = ml — im^tTt, in the top-quark propagator. At, resums leading finite-width effects and 
can be related order- by-order in as, aew to the renormalized mass and self-energy of the 
top quark. This is explicitly shown below in Section 13.31 The amplitude Ar^-^l is formally 
gauge independent, up to terms suppressed by 6 or higher. These gauge- violating terms, 
which would normally be cancelled by the electroweak background diagrams in Figure [11 
are numerically small as long as the condition 5 <^ 1 is satisfied. For the input parameters 
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and cuts adopted in Section HJ they amount to a correction to the leading resonant contri- 
bution to the cross section of much less than 1%. However, when S >1 the effective-theory 
expansion breaks down and these parametrically suppressed contributions can be numeri- 
cally important. In this case it is necessary to calculate the complete gauge- invariant set 
of diagrams, shown in Figure [H 

3.2 Virtual corrections 

The relevant one-loop contributions to the amplitude up to 0{6) are shown in Figure O 
As mentioned in Section I2.2[ we use the method of regions [35] to compute these loop 
integrals. This automatically yields the expansion of the full integral in 6 and has the 
advantage of simplifying the calculation. In the case of interest here, the two relevant 
momentum regions are soft (go ~ "f ~ ^) and hard (go ~ (f ~ !)• Only the hard part of the 
electroweak self-energy diagram (jJb) contributes to the amplitude at order 6, and this is 
automatically included in the cross section through the resummation of finite- width effects 
in the top-quark propagator, as explained in Section 13.31 Therefore, we focus here on the 
remaining QCD diagrams, (|2^) and ([2t)-(|2]F). 

To illustrate how the expansion by regions works we will now explicitly show how it is 
applied to the case of diagram (EJi), starting with the computation of the soft part. The 
full expression for the three-point loop integral is 



.(3,2) ^ ^ f d'^q ■■■i:^t-mth^{^t-(/-mt)--iv^2-(f)lMP2) /oox 

^(-1),. o^s^Fj ^^^y q2(^q2_2p^.q + D,)iq^-2p2-q) ' ^"^ ' 

where the ellipses denote constants and spinor quantities which do not depend on the loop 
momentum, q. Let us first consider the expansion of the three propagators appearing in the 
denominator. The gluon propagator has an homogeneous scaling, since all the components 
of q are of order S, and is not expanded. In the light-fermion propagator, g^ — 2p2 ■ g, the 
quadratic term (~ 6^) is suppressed compared to the linear term (~ 6) and can be dropped. 
Therefore, g^ — 2p2 ■ q — > —2p2 ■ q. In the top-quark propagator, the term quadratic in g 
can again be neglected, whereas p2 ■ q and Dt both scale as 6. Furthermore, given that 
Df — 2pt ■ g ~ (5 in the soft region, finite-width effects must be resummed to all orders in 
the top-quark propagator, leading to 

^ 1 . ^ 

^-. -n • (3-3) 



q^ -2pf q + Dt Aj - 2pt ■ q 

In the numerator of Eq. 03. 2p , the loop momentum, g, is always parametrically smaller than 
the external momenta, Pi, and can be neglected at leading order in 6. From the properties 
of the Dirac algebra and of the spinors, u{p), it also follows that 

^21MP2) = '^P2,MP2) , (3-4) 

Wt - mt)YWt - rrit) = 2p'^{^t - mt) - Da'' ~ 2p^t{y(t - rrit) , (3.5) 
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where we have used Dt ^ S to drop the second term. Thus, the leading soft contribution 
to the loop integral (13. 2 p is 



.{3,2)5 



asCF{4:p2 -pt) 



d'^q 



mt)...u{p2 



{27iYq\At~2pfq){-2p2-q)' 



(3.6) 



Note that the leading soft part in the expansion by regions is equivalent to the well known 
eikonal approximation, 



{^2-</)Y 



-)■ 



P2 



(3.7) 



{p2-qY ' {-P2-q) 

Eq. (13. 6p is much simpler than the original integral (13. 2 p and can be easily computed with 
standard techniques. The explicit result is given in Eq. (13. 8p . The parametric scaling of 
the correction (13.61) can actually be checked before the explicit calculation of the integral. 
Given that a soft momentum scales as 6; the gluon propagator scales as 5~^, the light-quark 
and top-quark propagator as 5~^ and the infinitesimal volume element d'^q as 6'^. Thus, 
Eq. (13. 8p is suppressed compared to the leading amplitude by as x 5^ x 5^"^ 5^^ 6^^ ~ 5^^"^, 
as expected. 

The remaining QCD diagrams in Figure [2] can be expanded in the soft region anal- 
ogously to diagram 1^). As a consequence of the simple structure of the quark-gluon 
vertex in the soft limit, Eq. (13. 7p . the soft corrections factorize in terms of scalar functions 
and the leading tree-level amplitude, ^[jij- The contributions of the different diagrams in 
Figure |2] read as follows (/i^*^ = e'^'^'^ fi'^'^ / {AtcY with /i being the renormalization scale): 



A 



(3,2)S 
(-l),a 



4(3,2)5 



—IGTriasCpfi 



r.2e"^? 



d'^q 1 



27r 



+ 2 



At J {27iYq^At-2pfq_ 

-2e 

^(3,0) 



A 



(3,0) 
(-1) 



A, 



jinit 



— IGniasCpfi "(pi ■ ps 







'(~i) 



d'^q 1 



(271)^^ g2 (_2p^ . q) (_2p3 . q) 



4(3,0) 



A 



(3,2)5 
i-l),d 



A 



(3,2)5 
(-l),e 



-IQiriasCpl^ \p2 ■ Pt) 



27r 



2^ + 24^^ 



OsCf 

2tt 



2^ + 24^^ 



d'^q 1 



1 



1 



Ai 



(271)"^ g2 (_2p2 .q)At- 2pt ■ q 

2e 



firrit 



A 



(3,0) 
(-1) 



-lQniasCFfi^''iP4 ■ Pt) I 



d'^q 1 



A, 
^,mt 



{2'kY q^ (-2p4 ■q)At- 2pt ■ q 

2e 



A 



(3,0) 
(-1) 



A 



(3,0) 
(-1) 



A 



(3,0) 
(-1) 
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.(3,2)5 
"^(-1),/ 



-IGniasCpfi "{pi ■ Pi 



cfiq 1 



1 



{2TiYq^ 



agCf 
2tt 

+Li2 






{S2t-mf){su-mf) 



2p2 ■ q) 

2\ 



-2p4 ■ g) At - 2pt ■ q 



A 



(3,0) 
(-1) 



mls2A 



m^S24, 
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-7r 



A, 



fimt 



-2e 



A 



(3,0) 
(-I)- 



(3i 



In the soft limit, diagram (2c) reduces to a scaleless integral that vanishes in dimensional 
regularization. This is a consequence of the fact that the upper quark line of the tree-level 
diagram ([T^) does not carry any information about the off-shellness of the intermediate 
top-quark propagator, and is thus not sensitive to the soft scale, mt6. The total soft 
correction is 



A 



(3,2)S 
(-1) 



|(3,2)S 



i 

a.^Cp ( Ai 



bV'Af^, 



271 



+ 2 



+ Li2 



-2e 



1-ln 



{S2t-m^t){su-mj) 



mls24 



{S2t-m^)is4t-mf) 



m1s24 



(3.9) 



In the hard region, the loop momentum scales as g ~ 1 and cannot be neglected in the 
numerator or denominator of Eq. (I3.2p . In this case, however, p^ — mf ~ 5 is much smaller 
than any other invariant and one can expand the integrand of Eq. (13. 2p in Df. The leading 
hard contribution for diagram ([2li) then reads 



A 



{3,2)H 
(-l),d 



agCp 



d'^q ...{^t - mt)-f^'{^t - 4- 'mt)...{^2 - (t)liju{p2) 



(2vr) 



g2(g2 _ 2p^ . g)(g2 _ 2p2 ■ q) 



(3.10) 



As for the soft region, the expansion in 6 leads to a significant simplification of the inte- 
grand. Note that, in the hard region, there is no resummation of self-energy insertions, 
since by assumption g^ — 2pt ■ q ~ rnf. Furthermore, at leading order in S, the hard part 
of (13. 2p coincides with the one-loop vertex correction to the on-shell single-top produc- 



tion process u{pi)b{p2) — )■ d{ps)t{pt), with p^ 



m 



f 



In the effective-theory language this 



is equivalent to the statement that the hard matching coefficients are obtained from the 
computation of fixed-order on-shell matrix elements. Strictly speaking, in this context the 
term "on-shell" would imply pf = fif. However, this condition must be satisfied only order- 
by-order in Og ~ Oew ~ S, and for an 0{6) calculation we can simply set pf = m^ in the 
leading hard contributions. Using the counting scheme (12.60 we can again determine the 
scaling behaviour of Eq. (13.101) . which is, in this case, a^ x 1 x 1 ~ 5^^'^, since all momenta 
scale as ~ 1. This is confirmed by the explicit result given below in Eq. (I3.12p . 
Applying the hard-region expansion to diagrams ([2]3)-([2]F) we obtain: 



A 



{3,2)H 
(-l),c 



^T/H .(3,0) 
^'^13^(-1) 
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.(3,2)// 



4(3,2)// 



sv.Ut'll + 



9e 
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(3,2)// 



(~i) ■ 2 

rT^// .(3,0) ^ 





vr [52](46)(3|2|l]a,CF 



m 



MwTw (si3 + M^)At 27r 2m2 - sst 
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vr 



[12](43)(6|4|5]a,Cj 






Mw^^w (si3 + M2,)At 27r 2ml - s^t 



In 



In 



S2t 



7711 



mt 



Sit - m^ 



mt 



'(-1),/ 
where the three scalar function 5V^, ^^it 1 ^^S ^^^ given by: 
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In 
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'^Kf 
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^^2? 



S2t 
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- mf 
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+ ^\n'{^] + 



S2t 



mt 



4(2m2 - 52*; 



S2t 



TTi: 



~2 



71 

24 



mtfi 



S2t 



ml 



2{2ml - S2t) 



In 



m. 



/i- 



(3.12) 



S2t— i'S4t 



The expressions in Eq. (I3.12p . once renormahzed as explained in the next section, agree with 
the results available in the literature (see e.g. Eq. (19) of Ref. |9] and Eq. (9) of Ref. [IH])- 
This is expected, since the leading hard contributions must coincide with the corresponding 
corrections to the on-shell top production process. The value of the parameter Xgc in 
Eq. (I3.12P depends on the particular version of dimensional regularization used for the 
calculation (see e.g. Ref. |51]), and is given by Xsc = in the 't Hooft-Veltman scheme (HV) 
and Xsc 



1 in the four-dimensional helicity scheme (FDH). As anticipated in Section [2^21 
the 0{as6^f'^) hard contribution of diagram ([2]f) vanishes since the intermediate top-quark 
propagator is off-shell by an amount ~ rn^. This leads to a further suppression in 5 
compared to the other one-loop QCD diagrams shown in Figure |2l In a strict effective- 
theory approach, the contribution of diagram (I2]F) would be encoded in an effective higher- 
dimensional hW^ — )■ hW^ vertex. 

The expansion of the self-energy diagram (JJ^i) in the hard region presents some sub- 
tleties. Let us consider the leading (fixed-order) top-quark propagator, i{j/t + mit)/Dt, and 
the QCD self-energy correction to the propagator. 



i{j/t + mt) 
Dt 



-ATiasCppi' 



r,2e 






Dt 



(3.13) 



By applying our counting to Eq. (I3.13p . one expects the leading hard contribution of the 
QCD self-energy correction to scale as UsS'"^ ~ 5"'^/^, i.e. it is enhanced by 5^^^"^ with 
respect to the leading-order propagator, which scales as 5~^ . Hence, in this case we have 
to push the expansion by regions one order further to include all the terms suppressed by 
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(5^/^ compared to the leading tree-level contribution. Thus, the expansion of Eq. fl3.13p in 
Df, up to the relevant order, is 



27T 



Xq 



2e 2 



TTi; 



/i^ 



2im 



tVI/t 



+ nit] 



D? 



+ 



irrit 



+ nit] 



D, 



(3.14) 



As anticipated, the first term in Eq. fl3.14p is "superleading" and in a generic renormal- 
ization scheme has to be resummed in the top-quark propagator. However, in the on-shell 
scheme adopted in this paper, the superleading bit is cancelled exactly by the top-quark 
mass counterterm, as shown in Section 13.31 

As pointed out at the beginning of this section, only the hard part of the electroweak 
self-energy diagram ([2]d) is relevant to our calculation since a soft contribution would scale 
as S^/'^Uew ~ f^^^^- In this case, the dominant hard correction to the propagator scales as 
otew^^"^ ~ ^^^ and is resummed, while higher-order terms are beyond the target accuracy 
pursued here and can be neglected. In principle, diagrams analogous to ([2)d), but with a 
photon or a Z-boson in the loop, should also be resummed. As we will see in the next 
section, in the on-shell scheme only the imaginary part of the one-loop two-point function 
is effectively resummed in the propagator, and diagram (^p) is therefore sufficient. 



3.3 Renormalization and resummation of finite-width effects 

We now discuss the renormalization of the one-loop amplitudes computed in Section 13.21 
In this paper we adopt the on-shell scheme (see e.g. Refs. [221IS3])- In this scheme, the 
wave-function counterterms for external massless particles vanish identically, i.e. 5Zi = 
for i = 1, ...,6. Furthermore, the electroweak coupling, Qew, is not renormalized at 0{as)- 
Thus, the only counterterm relevant for our calculation is the top-quark mass counterterm, 
Srut. This induces a renormalization of the top-quark propagator, which is 



+ mt) 



D, 



-idnit 



+ mt) 



Dt 



mt 






with the 0{as) mass counterterm given by 



+ 



zrrit 

'a 



(3.15) 



mt 



UsCb 

27r 



— + 2 + — 
2e 2 



m. 



fi^ 



(3.16) 



Comparing Eqs. ( 13.14^ and (]3.15p . it is clear that the superleading terms (and part of the 
subleading ones) arising from the self-energy correction are cancelled by the counterterm. 
Rewriting Eqs. (13.141) and (13.151) as contributions to the amplitude, one obtains 



j^i5,m ^ ^(3,2),rcn 



\-l),a 



'(-1) 



asCf 

2n 



^2^~ ~ ~2 



°^t \ 4(3,0) 



(3.17) 



In an arbitrary renormalization scheme, which in the following we will denote by r, there 
is, in general, no exact cancellation of the superleading contributions in Eqs. (I3.14p and 
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(I3.15p . and these terms are resummed inside the complex pole of the top-quark propagator, 
fif = mf — irritTt. Consider the renormalized top-quark two-point function, 

Y,(j/)=^-mt,r+%,M. (3.18) 

where rrit^r is the renormalized mass in the generic scheme r and J^t^r represents the renor- 
malized top-quark self-energy. The complex pole, /it, is defined as the solution for ]/ at 
which Eq. (13.211) vanishes (see e.g. Ref. [51]). Since At ~ 6, Eq. flS.lSp must include all 
contributions to Sj^^ up to 0{6), i.e. the one-loop QCD (~ as ~ S^^^) and electroweak 
(~ aew ~ S) top-quark self-energies, and the two-loop QCD contribution (~ a^ ~ 6). 
Other corrections to Sj ,. can be included perturbatively and need not to be resummed. 

In the following we will adopt the notation of Ref. [SS] and parametrize the top-quark 
self-energy as 

%M = ^IriP^yPL + ^lAP^yPR + KiP^)mt,r • (3.19) 

The renormalized quantities in Eq. (I3.19P are related to the unrenormalized ones by: 



tM) - 


- ^i:ip') + 5Zt, 




Kip') - 


- Sf(p^) + 5ZJ 




Kip') - 


= Ef(pV'"*^'' 


o 



(3.20) 
mt,r -l 

where bmt:!- and bZt^^ denote the mass and wave-function counterterms, in the scheme r, 
respectively. Eq. (13.181) can be rewritten in a form more suitable to the extraction of the 
pole by inserting identity matrices, in the form I = {]/ — mt^r){j/ + 'fnt,r) I {p^ ~'^'tr)i ^^ both 
sides of Tit,r and expanding the quantity {]/ + mt^r)^t,r{^ + 'nT't,r) around p^ = m^^,. This 
leads to 

rt(2/) = 2/- mt^r + ml^Vtt^rip^)— Y~ 

p"= — m^,, 

2st (/) + 2^t{p^) + HUp')] 7^^^ + - , (3.21) 



+'^L 



- {p^ — m 



t,r, 



where the ellipsis indicates terms that are suppressed beyond 0{6) and, thus, are irrelevant 
to the calculation presented here. For later convenience, in Eq. (I3.2ip we have introduced 
the function Qfj. = E^ ^ -|- S^ j, -|- 2S^ ,.. 

The complex mass, /ij, is obtained by solving the equation Tt(fit) = 0, and reads 



fJ't = rrit^r ^i't,r-("^t,r) + ^~^V("^t,» 



(3.22) 

Eq. (I3.22P can be further expanded in a^ and aew To this end, we introduce a decompo- 
sition of the self-energies analogous to Eq. (12. 7p 

^S _ ^5,(0,2) 2v^5.(0,4) ^5,(2,0) 
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;=^y,(0,2) 2^^.(0,4) ^L/R,{2,0) 



'■s^t,r 



(3.23) 



Note that the pure QCD part of the self-energy does not have an axial component, and 



■=fL,(0,n) -=^^,(0,71) _ -=^^,(0,71) 
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for n = 2, 4. Using Eqs. fl3.22p and fl3.23p one can easily obtain 



the following results for the pole mass and width: 
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(3.24) 



where we have used the information that only S^^ ' has a non- vanishing imaginary part, 
and all the functions are evaluated at p^ = m 



t,r- 



In the on-shell scheme (os), Eq. (13.241) assumes a particularly simple structure: 
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(3.25) 



The result (I3.25P follows from the particular form of the mass and field counterterms in 
this scheme l53l: 
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(3.26) 



Thus, at the level of accuracy we are interested in, the pole mass, rrit, can be identified 
with the on-shell mass, mt^osU- Furthermore, as anticipated at the end of Section [321 oiily 
the imaginary part of the one-loop electroweak self-energy is effectively resummed in the 
propagator. In particular, in Eq. (I3.25P the dependence from the one-loop (and two-loop) 
QCD contributions, which are responsible for the superleading terms in Eqs. (I3.14p and 
(I3.15p . cancels completely. The imaginary part of the complex pole is given, as expected, 
by the on-shell top decay width 
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(3.27) 



The residue of the propagator at the complex pole, 6Rt, which enters Eq. (12. Sp . can be 
easily extracted from Eq. (I3.2ip and expanded in as and aew up to the required accuracy. 



•^It is well known that the relation p.25p is corrected by terms of order a^^, which are, however, beyond 
the accuracy pursued here. 
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Here we only mention the fact that in a generic scheme SRt^r ~ «« + ..., and the term 
5Rt,T ^(li) is parametrically of the same order as the one-loop QCD corrections computed 
in Section 13.21 Thus, it also contributes to the amplitude at order 5. However, in the 
on-shell scheme, the expansion of SRt^os starts at order aew and contributes a correction to 
the amplitude which is beyond the accuracy pursued here. 

We would like to conclude this section with a remark on how a "bad" choice of renormal- 
ization scheme can lead to a breakdown of the effective-theory counting scheme, Eq. (12. 6p . 
Throughout this paper we assumed the scalings Dt = pi — mf^ ~ 6 and At = pf — fif ^ 6 
for the bare and resummed top-quark propagator. This is consistent in the on-shell scheme 
as n1 — rnf^g = fif — mf ~ 6. However, in a generic renormalization scheme the two con- 
ditions are incompatible since fj.f — m1^ ~ a^ ~ y5. This is a problem well known in 
the context of applying effective-theory methods to quarkonium physics. Thus, while in 
principle one could choose an arbitrary renormahzation scheme (e.g. MS), in practice this 
can lead to complications and loss of transparency in the expansion in 5. In other words, 
the effective-theory approach adopted here naturally identifies a class of "good" renor- 
malization schemes defined by the condition rrit^r — rrit ^ S, of which the on-shell scheme 
represents a particular example. 

3.4 Real corrections 

The last class of 0{5^^'^as) contributions to the amplitude we have to consider are rep- 
resented by real gluonic corrections to the tree-level process (12.11) . i.e. u{pi)h{p2) — )■ 
d{p3)b{pi)l~^ {p5)i>i{pQ)g{pi) . The relevant Feynman diagrams are shown in Figure [31 The 
(leading) real-emission amplitude can be written as 

^(-1) (97) = S2,ni9L9sA^^i^ {9^) + 6r,T:^9L9sA^{ii^ (9^) , (3.28) 

where gf denotes the two possible helicity states of the emitted gluon, and AL, and ^4^21 
represent the contribution of diagrams with a gluon attached to the upper or lower fermion 
line, respectively. The helicity amplitudes for the upper-line emission read 

,(3,1). +. _ r^^v^(46)(3|4 + 6|5](3|7-l|2] 

^\31] ydl ) — 9e 



[31] vyr ; ye.y ^^^^ ^^^^^ ^ ^2^^^^ ^^^^^3^^ 



vr 



V2 (46)[12][1|3 + 7|4 + 6|5] 



"^^^'"^^^'^ ~ ^HmwTw ^ "(.137 + M^)Ai [17][37]" ' ^^'^^^ 

where S137 = S13 -|- S17 — S37. For the lower-line emission we obtain 

(3|4 + 6|5](4|7-2|1] [12] (3|1 + 2|7](4|6|5] - /x,^(34)[57] \ 
(27)(47) ^(47) A„ ) 
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•(3|1 + 2|5](6|4 + 7|2] (46) (3|1|2](7|4 + 6|5] + M?(37)[25] 



[27] [47] [27] 



,(3.30) 



with At7 = {pt + ViY ~ lA- The complex mass, /zj, in the numerator of Eq. (I3.30p . which 
follows from the on-shell matching condition, p^ = fif, guarantees that QCD Ward identities 
are satisfied exactly. As already pointed out for the hard virtual corrections computed in 
Section [321 the matching condition needs to be satisfied only order-by-order in 6, and one 
could equally well replace fit with rrit in the numerator of Eq. fl3.30p . This would lead to a 
gauge-invariance violation proportional to Tt/rrit ~ S, which is a higher-order effect in our 
counting scheme. 

Besides the process ub — t- dbl^uig, at order 5^^'^ as the gluon-initiated processes 
ug — >• dbl^ Pib and gb — )■ dbl^ uiu also contribute to the amplitude for single-top pro- 
duction. The corresponding helicity amplitudes can easily be obtained from the crossing 
of Eqs. (13:291) and ^M)- 

3.5 Comparison of various approximations 

The results presented in this section have been obtained by making kinematic approxima- 
tions to the virtual and real corrections. Because real and virtual corrections are infrared 
divergent, we have to make sure that the approximations are consistent in that the real and 
virtual soft and coUinear singularities still cancel. Since, in Section HI we will be comparing 
our results to previous calculations, it is instructive to discuss the cancellation of infrared 
singularities and the relation between successive approximations to t-channel single-top 
production. 

The first approximation we consider is to treat the top as a stable particle [71jTmi26]. 
We will call this the stable-top calculation and denote the process hj ub ^ dt. Within 
this approximation, the renormalized virtual corrections read 

AZ\,t = AT,%,t (svi + 5Vi + 1^) + finite , (3.31) 

where A^^l^_^^^ is the corresponding tree- level amplitude; 5V^, 5V^ and 5mt/mt are given in 
Eqs. fl3.12p and fl3.16p : and "finite" represents non-singular terms that cannot be factorized 
in terms of the leading-order amplitude. Thus, the virtual corrections in Eq. f l3.3ip are 
precisely the leading hard corrections of diagrams (c) and (d) in Figure 12] without the top 
decay, plus the renormalization of the external top-quark line, 5Zt/2 = 5mt/{2mt). The 
corresponding infrared singularities are cancelled, in the usual way, by the real corrections 
due to the process ub ^ dtg. 

The next approximation we consider is the calculation presented in Ref. [IB], which 

we call the on-shell calculation and denote by m6 — )■ dbW . In this computation, the 
production of an on-shell top quark is combined with the decay of an on-shell top. This 
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leads to additional singularities in the virtual matrix element, which is given by 

ubAdbw ubAdbw \ ^-^ ^^ ^^ rrit J 

The additional contribution, 5V^ (given in Eq. fl3.12p ). corresponds to the one- loop cor- 
rection to the top decay, while 6mt/mt = 2 x 6mt/{2mt) accounts for the renormalization 
of the top-quark legs in the production and decay vertices. In order to match this with the 
real corrections, the latter have to be treated with some care. In particular, the interference 
between real gluon radiation from the production part of the process with the decay part 
has to be neglected. However, real corrections restricted to either the production or decay 
part are taken into account and the corresponding singularities cancel those of the virtual 
corrections, Eq. f l3.32p . 

In the calculation presented in this paper we move from on-shell top production to 

resonant top production, indicated by a star, ub ^ db W. As discussed in the previous 
section, this involves a whole tower of additional contributions. However, to the accuracy 
we are aiming for, the only additional contribution is due to the (leading) soft part of the 
diagrams in Figure [21 given in Eq. (13. 8p . Thus, we have additional virtual singularities, 
given in Eq. (13. 9p . and we obtain 

<-..w = ^'?*..w f ^^13 + SV,^ + 5Vi + ^ + dV") + finite , (3.33) 

with J^'^'^^ = A/^A. The new contribution, 5V^, exactly cancels the single poles in 
Eq. (I3.32P due to soft emission off a massive top-quark leg. Thus, Eq. (I3.33P contains 
only singularities corresponding to coUinear or soft emission off massless legs, as it should 
for a process with only massless external states. The additional poles in SV^ are to be 
cancelled by real singularities which essentially correspond to those neglected in the pro- 
cess m6 — )■ dbW, i.e. the interference between real radiation from the production and 
decay of the top. From a more formal point of view, one has to compute the real matrix 
element squared of the process ub — ?■ dbW g and expand it in 6. Because we only need 
the leading term in 6, it is sufficient to consider the diagrams of Figure |3l Taking the 
amplitude squared corresponding to these diagrams, integrating over the phase space and 
expanding the result in 6 leads to real singularities that precisely cancel those of the virtual 
corrections to the process ub ^ dbW . 

We stress once more that the calculation presented here is only meaningful if the top 
is nearly on-shell. If we are interested in the process ub -^ dbW with no constraint on 
the invariant mass of the final state particles, we have to compute the full virtual and real 
corrections, i.e. also take into account the diagrams that have been omitted in Figures[2]and 
El Needless to say, such a computation is considerably more involved, but the cancellation 
of real and virtual singularities works in a straightforward way. 
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4 Results 

The results for the three approximations outhned in the previous section, ub — )■ dt, 

t t* 

ub — )■ dbW and ub — )• dbW have been implemented using Monte Carlo integration, 

allowing us to calculate numerical values for both cross sections and kinematical distribu- 
tions. In this section we present a selection of these results and include comparisons to 
existing results where available. For simplicity, throughout this section the cross sections 
calculated using the three approximations will be referred to as o"P™^, a* and a** for the 
stable-top production, on-shell production followed by decay and resonant-top calculations 
respectively. 

4.1 Total cross sections 

We begin by presenting a comparison of our results, obtained using the three approxima- 
tions, to those for the production of a single, stable top quark as presented in Ref. [26] . 
We compare the results for an LHC run at a centre of mass energy y/s = 10 TeV, and 
use MSTW2008 PDFs [SS] and the corresponding strong coupling. Renormalization and 
factorization scales are set to be equal to a value of mt/2. The other input parameters 
used for this calculation are shown in Table [1] In the top-quark propagator, we use the 
tree-level decay width for the LO cross section and the ag-corrected width for the NLO cal- 
culation. This ensures that, at leading order, we obtain agreement between the stable-top 
production cross section and that of on-shell production plus decay, after integration over 
the fully inclusive decay of the top quark and IV-boson, (i.e. ctq'^'^'^-^ = rj"*-*). The results of 
our calculations using these parameters are shown in Table [2l We include the results for 
on-shell production taken from Table 5 of Ref. |26] for ease of comparison. 

The agreement with the existing results at both leading and next-to-leading order is 
very good when we perform the calculation for the production of a stable top. When we 
also include the subsequent decay of the on-shell top we still have good agreement at LO, 
but we no longer agree at NLO. This discrepancy is due to the use of the improved narrow 
width approximation for the on-shell top. We can see from 

prod decay ^prodpLO 

^0 = ^^LO = pLO = ^0 (4-1) 

that at LO the dependence on the decay width of the top quark cancels. However, at NLO 
there is no full cancellation and we are left with a residual dependence on the width. 



a 



prodpNLO _|_ ^prodpLO pLO _ pNLO 
t _ "O ^t ~rOl ^t _ _prod , prod^|_i ^ 



^prod ^ ar'^^-—^ . (4.2) 



pNLO '1 pNLO 



The difference between a* and a'^^°'^ is formally a higher-order correction, since a^™ (rj"*^ ~ 
r^^*-*) ~ a^, but leads to a visible numerical effect. From the resonant top calculation (last 
column in Table [2]), it is clear that taking into account the non-factorizable corrections 
causes a small, but noticeable increase of the cross section at both LO and NLO. 



'5 



23 



mt = 172 GeV 
Mw = 80.4 GeV 
aeyj = 0.03402 



LO 



r[ 



r 



NLO 

t 



2.05141 GeV 
1.46893 GeV 
1.32464 GeV 



Table 1 : Input parameters used for calculating the cross sections shown in Table [2j 



LO (pb) 
NLO (pb) 



Ref. 1261 



76.6 

84.4 



a 



prod 



a 



a 



76.62(1) 76.62(1) 77.36(5) 
84.41(1) 84.91(2) 86.3(3) 



Table 2: Comparison of total cross sections, calculated using our three methods, to those 
of Campbell et al. (2009) [26] at leading order (LO) and next-to-leading order (NLO). 



We now move on to discuss the comparison of our results to those in Ref. [18] for the 
on-shell production of a single top quark followed by its decay. We compare the results for 
an LHC run with centre of mass energy ^/s = 14 TeV. The input parameters can be found 
in Table IV of Ref. [18] . For this comparison we set the renormalization and factorization 
scales equal to mt and use the MRST2002 NLO PDFs [56] and the corresponding as value. 
The results are shown in Table [31 The total cross section in the stable-top production 
case is obtained by multiplying the production cross section, a^'^°'^, by the leading-order 
branching ratio for the top-quark decay, Bt^heu = 0.1104. 

As in the first comparison, we obtain good agreement between our results and the ex- 
isting results for the stable-top production and on-shell production plus decay calculations. 
The inclusion of the non-factorizable corrections again has the effect of increasing the cross 
sections by a modest amount at both LO and NLO. 

Finally, we look at the total cross section for our minimal realistic setup, described in 
Section 12.11 For illustration, jets are constructed using a standard k± cluster algorithm 
with the resolution parameter set to Drcs = 0.7, but any other jet definition could also be 
used. We assume that we can always identify Jb, the jet containing the b quark. We apply 



prod r) 



a 



prod rj 



t^fbeu 



On 



a 



On 



a 



ReL [IE] (pb) 
Our results (pb) 



17.69(1) 
17.71(11 



17.05(2) 
17.04(1) 



17.69(1) 
17.71(1) 



16.98(2) 
16.98(1) 



N/A 
17.94(1) 



N/A 
17.33(8 



Table 3: Comparison of total cross sections, calculated using our three methods, to those 
of Campbell et al. (2004) [18] at leading order (uo) and next-to-leading order (a). 
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mi = 172 GeV Mz = 91.2 GeV 

Mw = 80.4 GeV priJb) > 20 GeV 

ae«; = 0.03394 prie) > 25 GeV 

Tw = 2.14 GeV ^t > 25 GeV 

pNLO ^ 1,32813 GeV 120 < mj^y < 200 GeV 

Table 4: Input parameters used for calculating the cross sections shown in Table [5l 



LO (pb) 
NLO (pb) 






2.6786(1) 2.519(1) 
2.3079(1) 2.227(4) 



Table 5: Comparison of leading order (LO) and next-to-leading order (NLO) total cross 
sections for our minimal realistic setup. 



cuts on PT^Jb) and Pri^), the transverse momenta of J^ and the positron respectively, and 
on the transverse missing energy, ^t- In addition, we require that the invariant top mass, 
which is defined as 



m-y 



\/{p{Jb)+p{e)+p{u)f, (4.3) 



is close to the top-quark resonance. We perform the calculation for an LHC run with 
^/s = 7 TeV and use the MSTW2008 NLO PDFs. The renormalization and factorization 
scales are set to mj/2. A full list of the input parameters and cuts can be found in Table HI 
We use NLO PDFs and the NLO top decay width at all times to accentuate the off-shell 
effects by ensuring that the corrections are not coming from the change of PDF or width. 
As we now apply cuts to the decay products of the top quark, it is no longer possible to 
include the stable-top production calculation, as this requires fully inclusive decays. 

The results for the total cross section are shown in Table [51 With the introduction of 
the cuts, the effects of the non-factorizable corrections become more pronounced. Instead 
of increasing the cross sections, as was the case in the earlier comparisons (Tables [2] and 
[3]), the cross sections are now decreased by the inclusion of these corrections. The relative 
size of the off-shell contributions also increases from ~ 1% to a few percent. 

4.2 Distributions 

Partial results for kinematical distributions computed with the method presented here 
have been shown in Ref. [57] . In this section we will consider the minimal realistic setup 
described in the previous subsection, and use the parameters and cuts given in Table HI 
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Figure 4: Transverse top mass distribution for on-shell (dashed lines) and resonant (solid 
lines) top-quark production. LO results are shown in red, pNLO results (see text for expla- 
nation) in green and full NLO results in red. The orange line shows the subleading QCD 
contribution. 



The first distribution we present is that of the transverse mass of the top, defined as 



^T = \PT{Jb)\' + \pTie)\' + \pTiu)\^ - {pT{Jb)+PT{e)+M^)f ■ 



(4.4) 



We will compare on-shell {ub ^ db W) and resonant {ub ^ db W) top-quark production, 
shown as dashed and solid lines respectively in Figure HI To assess the importance of higher- 
order corrections, we will consider LO results (blue); full NLO results (red), including also 
partonic processes with gluons in the initial state; and partial NLO results (pNLO, green), 
including only one-loop corrections to partonic processes that are present at LO. Overall, 
the differences between the on-shell and resonant results are small. However, the off-shell 
effects do significantly change the shape of the distribution near the boundary m^ = rrit. As 
expected, the LO on-shell distribution shows a sharp edge at this boundary. The on-shell 
pNLO and NLO results have a contribution for r/ix > f^t, because the & jet, Jb, can contain 
gluon radiation, however, this contribution is very small. In the resonant calculation, we 
get a contribution in the region m-x > rnt even at LO. This contribution is reduced at NLO 
but is still significantly larger than in the on-shell result. 

It can be seen from Eq. (12. 8p that there are also terms of higher order in 6 in the squared 
tree-level matrix element. First, there are the interference terms gf^2Ile{A,2i) [^L) ]*)' 
second there are the terms proportional to g'^^g^lA^^ 



2)12 



QCD background; and, finally, there are subleading electroweak corrections g^„. \A) 



commonly referred to as the 

1(3,0) 1 2 
\0) I • 

The former two are 0(5^), whereas the latter is 0{6^). Therefore, parametrically, these 
terms are beyond our NLO approximation. However, they are a subset of subleading 
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Figure 5: -ffT(had) distribution for on-shell (dashed) and resonant (solid) top-quark produc- 
tion at LO (blue) and NLO (red). The orange line shows the subleading QCD contribution. 



corrections that are very easy to compute and it is useful to compare their numerical 
importance with respect to the terms of 0{5^^'^). Thus, the QCD background is shown 
separately in Figure HI This contribution is actually important in the region m^ <C m^, 
but is insignificant when mx ~ rnt- The subleading electroweak terms of 0{5'^) and 0{5^) 
are very small, and on the scale of Figure H] they would show simply as a straight line along 
the bottom of the plot. 

We now move on to look at our second distribution, the sum of the hadronic transverse 
momenta, defined as 

HT:ih^d) = \pT{Jb)\ + \pT{J{)\ , (4.5) 

where Ji is the (non h) jet with the largest transverse momentum. The results are shown 
in Figure [51 In this case, the corrections due to the off-shellness of the top are relatively 
large. The on-shell and resonant results differ by up to 10% in some bins. On average, the 
difference is approximately 3-4%, in line with the results presented in Table |5l We have 
studied several similar distributions and the effect of the off-shell corrections is typically 
somewhat smaller than for i7T(had). Again, the QCD background is separately shown in 
Figure [51 As for the transverse mass, it is important only at the edge of the distribution, 
and the subleading electroweak effects are again too small to be shown. 

Finally, we turn to the invariant mass of the top, as defined in Eq. (14. 3p . The upper 
panel of Figure [HI shows the result at leading order for our resonant calculation, along 
with the next-to-leading order results for both the on-shell and resonant calculations. It 
should be noted that the leading order on-shell production and decay distribution is a delta 
function centred at the top mass and so would appear only as a single point at rrit in the 
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Figure 6: Upper panel: Invariant mass distribution for LO (blue) and NLO (red) resonant 
top production, compared to NLO (red, dashed) on-shell top production. Lower panel: 
Ratios of various corrections (full NLO correction (red), QCD background (orange), off-shell 
corrections (green), subleading electroweak corrections (azure)) to the full NLO result. 



figure. Therefore, it has been omitted. 

We can see from the figure that the inclusion of NLO effects causes a deviation from the 
Breit-Wigner shape of the distribution. We also note that there is a visible difference be- 
tween the NLO distributions obtained via the two methods. This difference is particularly 
noticeable at invariant mass values greater than rrit. 

Considering the shape of the invariant mass distribution also helps us to understand 
why the inclusion of off-shell corrections increases the total cross section when no cuts are 
applied but decreases it when we apply some cuts. At LO, the invariant mass distribution 
has a delta-spike shape in the on-shell case but a Breit-Wigner shape in the resonant case. 
Taking into account NLO corrections, this picture is modified somewhat, but comparing 
the on-shell and resonant distributions, the former still has a more prominent peak at 
minv = ^TT't, whereas the latter is larger in the region minv > '^t- Thus, there are two 
competing effects. If the cut on minv is mild enough, the increased contribution of the 
resonant result for rriinv > frit outweighs the larger peak of the on-shell result, whereas for 
a tight cut on minv, the on-shell result is larger. 

Note, however, that we must apply a cut to the invariant mass, otherwise our assump- 
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tion that the top quark is close to resonance is no longer valid and the (5-counting no longer 
applies. As mentioned before, if we take our power counting at face value we would expect 
our approximation to work in a window of rrimv — iTit ^ Tt- However, the effective theory 
actually works in a considerably larger window. To illustrate this, let us consider the lower 
panel of Figure [6], where we compare the relative importance of the various corrections to 
the full NLO resonant distribution. The NLO corrections, (Tnlo (shown in red), including 
all corrections of 0{S^^'^), are dominant in the vicinity of wiinv — rrit, as expected. They 
are larger than the (9(5^) QCD corrections, ctqcd (shown in orange), and the subleading 
electroweak corrections, ctew (shown in azure). However, away from resonance the for- 
mally subleading QCD corrections actually become numerically more important than the 
formally leading corrections. This is a clear indication that our power counting is no longer 
applicable in this region and, therefore, the effective theory breaks down. A similar point 
can be made in the case of the subleading electroweak corrections, which have been di- 
vided into (!?((5^) contributions (shown as solid azure line) and 0{6^) contributions (shown 
as dashed azure line). These corrections are much smaller than our NLO corrections. In 
the resonance region, the 0{6^) corrections are larger than the 0{6^) corrections, but for 
m-inv ^ 130 GeV this is no longer true. Once more, this indicates the limitations of our 
power counting in this region. We should also mention that for rrii^v ^ 160 GeV the NLO 
corrections are huge compared to the LO result and, therefore, our result is not reliable. 
Finally, the off-shell effects, defined as the difference between the full NLO results for 
the resonant and on-shell calculations, are shown in green in the lower panel of Figure El 
These effects are relevant near and above the resonance region, but are very small below 
resonance. 

We stress that the distributions presented here are only a sample of the types of dis- 
tribution that could be calculated. In principle, any infrared safe quantity with arbitrary 
cuts on the final state particles and jets could be easily computed. 

5 Conclusion and outlook 

In this work we have presented a method which allows the inclusion of off-shell effects in 
resonant-particle production with a minimal amount of computation. The method is based 
on a simultaneous expansion of the cross section in the couplings and the small kinematic 
variable At/rrit ~ Tt/mt. It has been applied to t-channel single-top production at the 
LHC The calculation includes the first non-trivial corrections to the narrow-width ap- 
proximation, corresponding to production-decay interference terms, and generalizes earlier 
results of one-loop corrections to t-channel single-top production. 

Generally speaking, off-shell effects are small for inclusive quantities. For the total 
cross section, for example, we find an effect of the order of 1%. However, depending on the 
cuts applied, off-shell effects can be sizeable. For most distributions we have considered, 
the off-shell effects amount to a few percent of the LO result, reaching up to 10% in more 
extreme cases, such as the H^ihad) distribution defined in Eq. (14.51) . In particular, they 
can significantly change the shape of distributions near phase-space boundaries related to 
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off-shell effects, such as the edge for the transverse mass. This is, of course, not surprising 
since sharp edges in distributions are usually related to having particles on-shell. Thus, 
allowing the top quark to become slightly off-shell can have a large impact in this region. 
The calculation presented here can also be seen as a proof of the workability of the 
effective-theory method. As we have shown in this paper, the inclusion of the leading off- 
shell effects is relatively straightforward, requiring in principle only the calculation of simple 
soft corrections, since the hard part of loop integrals can be easily related to results for on- 
shell production and decay of the massive particle. Furthermore, the computation of real 
corrections requires only minor modifications to the standard subtraction procedure. Given 
the generality of the effective-theory approach, the method can be easily applied to other 
processes of phenomenological interest at the LHC One such process is clearly represented 
by top-quark pair production, which will be intensively exploited for measurements of the 
top-quark properties, and whose study we reserve for future publications. 
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